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Closely  spaced  thin-wire  antennas  are  analyzed  by  the  moment  method 
using  the  piecewise-sinusoidal  function  to  describe  the  current  variation 
along  the  length  of  the  wire  and  a Fourier  series  for  the  circumferential 
variation.  Data  are  presented  for  the  circumferential  variation  of  the 
surface  current  density  of  parallel  dipoles  obtained  by  the  'thin-wire* 
theory,  wire-grid  model,  and  the  present  formulation.  < 
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I.  INTRODUCTION 


A fundamental  approximation  made  in  the  "thin-wire"  theory  is  that 
the  surface  current  density  is  uniform  around  the  circumference  of  the 
wires.  In  order  to  make  this  approximation,  two  conditions  are  imposed 
on  the  wires.  First,  the  wire  radius  must  be  much  less  than  a free 
space  wavelength,  and  second,  no  two  wires  pass  within  a few  wire 
diameters  of  each  other. 

This  report  considers  the  problem  of  closely  spaced  three-dimensional 
electrically  thin  wires.  Smith  [1]  and  Olaofe  [2]  had  considered  the 
problem  of  closely  spaced  infinite  cylinders,  and  Smith  included  a 
modification  for  finite  wires.  The  solution  presented  here  is  a modi- 
fication of  the  piecewise  sinusoidal  reaction  formulation  for  "thin- 
wire"  structures  [3].  A Fourier  series  is  used  to  represent  the 
circumferential  variation  of  the  axial  component  of  the  current. 

Numerical  data  illustrate  the  effects  of  wire  radius  and  separation  on 
the  circumferential  variation  of  the  current.  The  data  are  compared 
with  the  results  of  "thin-wire"  theory  and  also  with  a wire-grid  model 
of  the  closely  spaced  wires. 


II.  THEORY 

The  reaction  concept  was  introduced  by  Rumsey  [4],  and  applied  by 
Richmond  [3]  to  derive  the  reaction  integral  equation  (RIE)  for  radiation 
and  scattering  from  thin-wires  and  conducting  surfaces. 

Consider  the  problem  of  scattering  by  a wire  structure  enclosed  by 
the  surface  S,  and  with  the  impressed  electric  and  magnetic  currents 
(^i,Mi)  confined  in  the  volume  V].  Let  the  ambient  medium  be  free 

space.  The  ej“^  time  dependence  of  all  sources  and  fields  will  be  sup- 

pressed. In  the  presence  of  the  wire,  (Ji,Mi)  generate  the  field  (£,ii). 
From  the  surface-equivalence  theorem  of  ^^chelkunoff  [5],  the  wire  can  be 
replaced  by  free  space  if  the  following  surface-current  densities  are 
introduced  on  the  surface  S 

^ = n X H[  (1) 

= E X n (2) 

where  n is  the  outward  directed  unit  vector  on  S.  By  defining  (^,Ms) 
as  in  Equations  (1)  and  (2),  the  total  field  inside  the  wire  is  zero. 

The  RIE  for  the  unknown  currents 


1 


(3) 


II  + n|  = 0 


where  (E  ,h"^)  is  the  field  of  test  sources  (^,^)  located  interior  to 
the  surTace  S and  radiating  in  the  homogeneous  medium.  Below  we  will 
consider  perfectly-conducting  wires,  and  thus  ^ = 0. 

Let  each  segment  of  the  wire  structure  have  a circular  cross 
section  and  at  each  point  on  the  surface  define  a fight-handed  orthogonal 
coordinate  system  with  unit  vectors  (n,$,i)  where  n is  the  outward 
normal  vector,  ^ is  directed  along  the  wire  axis  and 

♦ = A X n.  (4) 

At  this  point.  Equation  (3)  is  usually  simplified  by  making  the 
following  thin-wire  approximations: 

1.  wire  radius  a « A 

2.  neglect  integrations  over  the  end  surfaces  of  the  wire 

3.  the  ^ component  of  ^ is  zero 

4.  the  A component  of  ^ is  independent  of 

We  will  use  all  of  the  above  approximations  except  for  #4.  The  surface 
current  density  on  the  wire  can  then  be  written  as 

^(A,4,)  = A J^(A,4,)  . (5) 

Inserting  Equation  (5)  into  Equation  (3)  yields 


L 2v 

J^(A,^)  {^E'")ad*  dA 
J 5 

0 0 


J|  - M.-H'")dv 


(6) 


where  L denotes  the  overall  wire  length. 
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Equation  (6)  will  be  solved  via  the  method  of  moments  [6].  To  do 
this  is  expanded  in  terms  of  a finite  series  as  follows: 


(7) 


= L Jn  G 

^ n=l  " 


where  the  Jp  are  unknown  coefficients  and  the  Gn(A,4i)  are  a known  basis 
set.  Enforcing  Equation  (6)  for  N distinct  test  sources  yields  the 
following  system  of  simultaneous  linear  equations 


I J-  z 

n=l 


= V 
n mn  m 


n = 1,2,»»*N 


(8) 


where 


' I 

i i 


mn 


JJ 


6^(4..^)  (^e"’)  ds 


(9a) 


\ = jjj  (Ji’i"'  - Mi -if)  dv  (9b) 


and  where  the  integration  in  Equation  (9a)  is  over  the  surface  of  the 
n-th  expansion  mode.  The  accuracy  and  computational  efficiency  of 
the  solution  is  dependent  upon  the  choices  for  the  expansion  and 
testing  functions. 


In  defining  expansion  functions  we  choose  the  piecewise-sinusoidal 
function  [3]  to  describe  the  l variation  and  a Fourier  series  for  the 
(|i  variation.  The  expansion  is 


Thus,  the  Gn(A,(|>)  are  of  the  form  Fn(^),  Fn(t)cos'|>.  Fn(^)sin(|i,  Fn(A)cos2(l> , 
Fn(‘^ )sin2(Ji , •••  Fn(t)cosK((>  and  F^(t)sinK<t>  where  the  Fn(t)  are  the 
piecewise-sinusoidal  V-dipoles.  If  P V-dipoles  are  used  to  describe 
the  A variation  and  K terms  in  the  Fourier  series  are  retained,  then 
the  total  number  of  unknowns  will  be 

N = P{2K+1).  (11) 


The  testing  functions  were  chosen  identical  to  the  expansion  modes. 
Thus,  ^ = 0 and  ^ This  is  an  application  of  Galerkin's  method 

and  the  impedance  matrix  will  be  symmetric.  Our  initial  choice  for  the 
test  modes  was  filamentary  sources  on  the  wire  axis.  This  choice  led 
to  serious  problems  of  relative  convergence  [7]  and  was  abandoned  for 
the  true  Galerkin  solution. 

The  above  choice  of  expansion  and  testing  functions  has  several 
advantages.  First,  the  expansion  modes  are  placed  in  an  overlapping 
array  on  the  wire  so  that  continuity  of  current  is  enforced.  Evaluation 
of  the  elements  in  the  impedance  matrix.  Equation  (9),  requires  a 
quadruple  integration,  i.e.,  two  to  find  and  two  to  integrate  over 
the  surface  of  the  n-th  expansion  mode.  However,  since  we  employ  the 
piecewise-sinusoidal  functions,  only  two  of  these  integrations  are 
done  numerically.  Advantage  can  also  be  taken  of  the  orthogonality 
of  the  Fourier  modes  in  evaluating  terms  where  the  expansion  and  test 
modes  overlap. 


III.  NUMERICAL  RESULTS 

In  this  section  numerical  results  will  be  presented  for  the 
surface  current  density  of  parallel  center-fed  dipoles,  as  illustrated 
in  Figure  1.  Results  obtained  by  the  method  of  the  previous  section 
will  be  compared  to  those  of  "thin-wire"  theory  [8],  and  also  to  the 
results  of  a wire-grid  model  of  the  parallel  dipoles  [9].  In  the  wire- 
grid  model  each  thin-wire  dipole  is  modeled  by  about  sixteen  extremely 
thin  wires  uniformly  spaced  around  its  circumference.  The  dipoles  are 
center-fed  by  unit  voltage  delta-gap  generators.  Although  ^ (*■,<(>)  varies 
along  the  A coordinate,  only  ^ around  the  center  of  the  wire  will  be 
presented.  The  following  calculations  demonstrate  the  ability  of  the 
preceding  formulation  and  computer  programs  to  account  for  the  (^ 
variation  of  ^ for  closely  spaced  thin  and  moderately  thick  parallel 
dipoles.  The  surface  current  densities  are  .normalized  with  respect  to 
an  isolated  dipole.  In  all  the  computations  to  follow,  3 modes  were 
used  to  describe  the  longitudinal  distribution  on  each  X/2  dipole. 


4 


rT 


1 


I 


The  normalized  magnitude  and  phase  of  the  surface  current  density 
of  two  parallel  X/2  dipoles  for  various  wire  spacings,  c/a,  and  wire 
radii  of  lO'^A  and  10"^A  are  shown  in  Figures  2 and  3,  respectively. 
Referring  to  Equation  (10),  the  coefficients  of  the  expansion  modes  for 
the  above  cases  are  presented  in  Tables  1-A  and  1-B.  Note  that  since 
^ is  an  even  function  in  <(>,  all  sink<|>  terms  vanish.  Following  the 
notation  of  Equation  (10),  K will  denote  the  number  of  dependent  terms 
in  the  surface  current  expansion.  For  example,  K=0  implies  that  only 
the  constant  ("thin-wire"  approximation)  term  is  used,  K=1  implies  that 
the  constant,  cos<|>  and  sin(ji  terms  are  used,  etc.  In  Figures  2 and  3 note 
that  the  results  of  the  previous  section  are  in  excellent  agreement  with 
the  wire-grid  model  of  the  parallel  dipoles.  In  Tables  1-A  and  1-B  note 
that  the  coefficients  of  the  constant  term  for  the  K=0  case  is  very 
close  to  that  of  the  K=1  and  K=2  cases,  except  in  the  extreme  case 
where  a = 10"2x  and  c/a  =1.1. 

Figure  4 and  Figure  5 show  the  normalized  magnitude  and  phase  of 
the  surface  current  distribution  for  three  parallel  X/2  dipoles  for 
various  wire  spacings  and  wire  radii  of  10~^A  and  10"^X,  respectively. 

The  coefficient  of  the  expansion  modes  are  tabulated  in  Tables  2-A  and 
2-B  for  the  above  cases.  In  the  cases  where  the  wires  are  very  close 
together  (c/a=l.l  and  c/a=2.0),  the  coefficient  predicted  by  the 
case  K=0  is  in  error  by  20%  to  more  than  100%. 
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X WIRE  GRID  MODEL 


(d)  ya  = 50.0 

Figure  2.  (Continued) 


itinu 


TABLE  1-B,  Two  Parallel  \/2  Dipoles  Each  of  Radius  lO'^X 


Note  that  sink<{i  terms  vanish 


c/a 

expansion  mode 
used  (K) 

Constant 

COSiti 

cos24. 

0 

7.96/-15.7 

- 

- 

1.1 

1 

8.19/-13.7 

20.64/149.0 

- 

2 

8.19/-13.6 

5.10/157.4 

0 

CM 

1 

O 

CO 

- 

- 

2.0 

1 

8.32/-9.9 

21.47/145.5 

- 

2 



8.32/-9.9 

4.89/153.7 

20.0 

0 

18.38/7.0 

- 

. 

1 

18.35/6.7 

24.51/83.1 

- 

2 

18.35/6.7 

1.93/133.6 

0 

13.49/-26.9 

- 

- 

50.0 

1 

13.56/-26.4 

10.90/-146.6 

2 

13.56/-26.4 

0.82/-76.1 



wire  #U1.1 


wire  #2<1.1 


wire  #K2.0 


wire  #2<2.0 


TABLE  2-A.  Three  Parallel  A/2  Dipoles  Each  of  Radius  10  A 

Coefficients  of  the  Expansion  Modes  for  0$  at  the 
Center  of  the  Dipoles  #1  and  #2  (mA/meter) 

Note  that  sink4>  terms  vanish 


expansion  mode 
used  (K) 


Constant 


139.36/-29.2 


79.47/-27.0 


83.79/-27.1 


84.25/.27.2 


107.35/145.5 


13.26/-11.8 


5.73/17.0 


5.25/25.2 


101.48/-28.1 


85.41/-27.0 


85.78/-27.0 


85.78/-27.0 


33.54/1 36.0 
6.39/63.0 
6.53/70.5 
6.53/70.5 


180.71/149.0 


170.37/147.9 


168.61/147.9 


168.61/147.9 


cos2{fi 

cos34i 

21.18/153.5 

- 

20.36/154.0 

6.81/-38.4 

66.09/151.4 


68.60/151.1 


17.16/151.9 

17.16/151.9  1.22/166.0 


60.49/149.2 
60.46/149.2  0.0 
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TABLE  2-B.  Three  Parallel  A/2  Dipoles  Each  of  Radius  10"^A 


Note  that  sinki|i  terms  vanish 


The  results  shown  in  Figures  2-5  and  in  Tables  1 and  2 suggest 
opposite  conclusions  as  to  the  accuracy  of  the  thin-wire  approximation 
that  ^ is  independent  of  First,  Figures  2-5  clearly  show  that  Js 
is  strongly  <t>  dependent,  even  for  c/a  as  large  as  50  and  a as  small  as 
A/1000.  Thus,  one  is  led  to  the  conclusion  that  the  thin-wire  approxi- 
mation is  not  valid.  Next,  Tables  1 and  2 show  that,  except  for 
extremely  small  c/a  where  the  thin-wire  approximation  is  known  to  fail, 
the  independent  term  predicted  from  "thin-wire"  theory  (K=0  case)  is 
essentially  identical  to  the  (|)  independent  term  from  the  K=1  or  K=2 
cases.  Since  for  wires  whose  radius  a <<A  it  is  the  ip  independent  term 
which  dominates  the  far-zone  radiated  fields,  "thin-wire"  theory  will 
predict  the  proper  far-zone  fields.  Also,  if  the  device  used  to  measure 
input  impedance  is  sensitive  to  the  average  value  of  voltage/Jg  around 
the  circumference  of  the  wire,  then  it  will  also  depend  only  on  the 
(ji  independent  modes.  Thus,  one  can  be  led  to  the  conclusion  that  the 
"thin-wire"approximation  is  in  fact  valid.  This  discrepancy  can  be 
resolved  by  stating  that  the  "thin-wire"  approximation  is  not  valid  for 
predicting  the  circumferential  distribution  of  current  on  a general 
wire  structure,  however,  it  is  useful  for  computing  such  quantities  as 
impedance  and  far-zone  fields. 

At  this  point  it  is  reasonable  to  ask  the  question  "why  go  to  all 
the  time  and  trouble  to  determine  the  i dependence  of  ^ since  it 
ordinarily  does  not  affect  either  the  input  impedance  or  far-zone  fields?" 
The  answer  is  two-fold.  First,  "thin-wire"  computer  programs  are  in 
such  widespread  use  that  any  information  relative  to  the  basic  approxi- 
mations made  in  these  codes  is  of  value.  Secondly,  there  are 
quantities  which  the  p dependence  of  ^ do  affect.  For  example,  the 
conductor  loss  resistance  is  a strong  function  of  the  circumferential 
distribution  of  current  on  the  wire  [1].  A knowledge  as  to  how  this 
loss  resistance  varies  with  wire  radius  and  wire  spacing  is  important 
in  determining  the  efficiency  of  such  antennas  as  the  electrically 
small  multiturn  loop  antenna. 


IV.  CONCLUSION 

This  report  has  considered  the  problem  of  determining  the  circum- 
ferential current  distribution  on  closely  spaced  electrically  thin-wires. 
The  solution  is  a modification  of  the  piecewise  sinusoidal  reaction 
formulation  for  thin-wire  structures.  A Fourier  series  is  used  to 
represent  the  P dependence  of  the  wire  surface  current  density. 

Numerical  data  presented  illustrate  that  there  can  be  a considerable 
circumferential  variation  even  when  the  wire  to  wire  separation  exceeds 
several  wire  diameters.  This  can  be  important  in  determining  antenna 
efficiency;  however,  if  these  wires  were  to  be  analyzed  by  "thin-wire" 
theory  (i.e.,  neglecting  the  circumferential  variation  of  the  surface 
current  density)  the  correct  impedance  and  far-zone  fields  are  usually 
obtained. 
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